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In-plane angular dependent resistivity ADR was measured in the non-Fermi liquid regime of 
CeCoIns single crystals at temperatures T < 20 K and in magnetic fields H up to 14 T. Two scaling 
behaviors were identified in the low-field region where resistivity shows a linear T dependence, 
separated by a critical angle 6 C which is determined by the anisotropy of CeCoIns; i.e., ADR depends 
only on the perpendicular (parallel) field component below (above) 8 C . These scaling behaviors and 
other salient features of ADR are consistent with d-wave density waves. 



The recently discovered heavy fermion compound 
CeCoIngi has generated a lot of interest, partly due 
to the many analogies present between this compound 
and the high transition temperature T c superconduc- 
tors. Like the cuprates, CeCoIns has a layered 
crystal structure^, a quasi-two-dimensional electronic 
spectrum^, and a superconducting phase appearing at 
the border of the antiferromagnetic phased. NMR Knight 
shift** and heat conduct ivitj* measurements revealed 
even parity pairing consistent with d-wave pairing sym- 
metry, while specific heal*^ and electrical transport 
experiments have shown typical non-Fcrmi liquid behav- 
ior and have evoked the possibility of a pseudogap. In 
fact, infrared spectroscopy studies have revealed the de- 
velopment of a gap in the density of electronic states of 
CeCoIns below ~ 100 K 11 . As in the case of cuprates, 
all these features suggest the existence of quantum crit- 
ical phenomena 12 . However, the Fermi surface topology 
is much more complex than in cuprates, with multiple 
sheets, which implies the participation of several bands 
in the pairing processi^iii^. 

In this complex picture, the nature of the normal state 
in both systems is not trivial. Until recently, the nature 
of the pseudogap (or the non- Fermi liquid) phase present 
in CeCoIns has not been addressed. The giant Nernst 
effect found in this material above Tj&, however, re- 
opened the problem of the origin of the pseudogap state, 
not only in this compound, but also in the cuprates. 
It has recently been shown that the giant Nernst ef- 
fect in CeCoIns is consistent with unconventional den- 
sity waves UDW or d-wave density waves d-DW—. We 
recall that in underdoped cuprates, the pseudogap phase 
has been attributed to the existence of phase fluctua- 
tions of the superconducting order paramete:»i2ii2i 2 £i 2 i2 2 i 
as well as d-DWSiSiS. In particular, the giant Nernst 
effect observed in I^-zSr^CuO^, YBa2Cu307_a, and 
Bi2Sr2CaCu208+,5 has recently been interpreted satisfac- 
torily in terms of d-DW—. 

Here, we present new information concerning the na- 
ture of the pseudogap of CeCoIns through angular de- 
pendent resistivity ADR. Two different scaling behaviors 
are revealed in this compound in the non-Fermi liquid 



regime for applied magnetic fields less than 5 T. Specif- 
ically, for angles below (above) a certain critical angle 
6 C , the field component perpendicular (parallel) to the 
a6-planes determines the angular-dependent dissipation. 
The critical angle 9 C is determined by the anisotropy of 
the material. These results and the other salient fea- 
tures of ADR measured in high applied magnetic fields 
are consistent with d-wave density waves with Landau 
quantization of the quasiparticle spectrum. Our findings 
bring further understanding of the underlying physics in 
the non- Fermi liquid regime of this exotic compound and 
could also advance the understanding of the pseudogap 
state of the cuprates. 

High quality single crystals of CeCoIns were grown us- 
ing the flux method. Typical sizes of the crystals are 

0. 5x0. 5x0. 1mm 3 with the c-axis oriented along the small- 
est dimension. We determined the out-of-plane p c and 
in-plane p a b resistivities using the electrical contact con- 
figuration of the flux transformer geometry, as described 
in Ref;— . The angular dependent resistivity p a b{9) was 
measured by rotating the single crystal from ff | \c— axis to 
ff||f ||a-axis, with the angle between ff and the c-axis. 
Here we report only in-plane resistivity results. 

The angular dependent in-plane resistivity of CeCoIns 
single crystals was measured in the non-Fermi liquid 
regime, for T < 20 K and in magnetic fields ff up to 
14 T. Figure 1 is a typical contour plot of resistivity ob- 
tained from p vs 9 data measured at different fixed ap- 
plied magnetic field values and at a temperature of 6 K. 
A similar topologic contour is obtained over the whole 
T range investigated. The squares, circles, and triangles 
correspond to n = 0.45, 1, and 2/3, respectivelywhich 
consistent with previously report^ 2 *, where n is the ex- 
ponent of the power-law T dependence of the resistivity; 

1. e., p(T) — po + AT n . Note that the contour plot can be 
divided into three different regions based on the shape 
of the topologic plots and the values of n. Region I is 
the rectangular part, for which n — 1. This region is 
the so-called low-field region (ff < 5 T) for reasons given 
later on. Region II is the elliptical part, for which n 
= 0.45 and which corresponds to the region around the 
maximum value of the p{H) curves for different angles. 
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Region III is a large area with n — 2/3, which is outside 
of regions I and II at even higher fields. Region I is more 
interesting since it represents the typical non-Fermi liq- 
uid behavior (for this region n = 1) and displays the two 
scaling behaviors. Therefore, we will discuss the anoma- 
lous magnetoresistivity behavior in this region in more 
detail. 

The bottom Inset to Fig. 2(a) shows p(9) measured 
in a low magnetic field of 2 T. [Data were measured in 
region I.] The resistivity decreases as 9 increases from 
0°, reaches a minimum around 60°, and then displays a 
peak at 90°. The data of the top Inset to Fig. 2(a) were 
measured in a high field of 9 T. [Data were measured in 
region III.] The resistivity increases with increasing angle 
and displays a shoulder, followed by a peak at 90°. 

In the low field region (region I), a spectacular scaling 
of these data is achieved by plotting the resistivity as a 
function of the component of the applied magnetic field 
perpendicular to the afo-planes, i.e., H± = HcosQ, as 
shown in the main panel of Fig. 2(a). Note that all 
the curves measured in H < 5 T overlap for all 9 values 
between 0° and a critical angle 6 C w 60°, marked by the 
arrows, and deviate from this H cos 9 scaling at higher 
angles. This value of the critical angle is independent 
of temperature and applied magnetic field for the range 
investigated. 

We employed a second protocol to measure resistivity, 
in which we kept the angle between the magnetic field 
and the c-axis fixed and scanned the magnetic field. Fig- 
ure 2(b) shows the resistivity as a function of the field 
component parallel to the a— axis (H\\ = Hsin9), mea- 
sured at 6 K for different orientations of H. Note that 
all the resistivity curves scale this time with Hsin9 for 
9 > 9 C ~ 60°. Therefore, the resistivity data follow two 
scaling laws: 



P (H,9) 



fi(Hcos9) for 9 < <9 C « 60° 
f 2 (Hsin9) for 9 > 9 C w 60°. 



(1) 



Figure 1 also clearly shows the presence of the two scaling 
laws in region I as evidenced by the rectangular region of 
the p contours. The presence of these two scaling behav- 
iors indicates that below (above) the critical angle, the 
field component perpendicular (parallel) to the a6-planes 
determines the in-plane dissipation. A similar analysis of 
the angular dependent resistivity data measured in high 
magnetic fields (region III) has shown that the data are 
still dominated by the perpendicular (parallel) field com- 
ponent below (above) the critical angle, although the two 
scaling laws are no longer present. 

Next we try to understand the above shape of p(9) 
measured in both low and high fields and the relation- 
ship between the two scaling behaviors. We plot in the 
Inset to Fig. 3 the resistivity measured at 6 K in scan- 
ning H up to 14 T at the two fixed angles (9 = 0° and 
9 = 90°) corresponding to the two field orientations (per- 
pendicular and parallel, respectively, to the a6-planes) 
which seem to determine the physics below and above 
9 C , respectively. For the first field orientation (H _L ab 



planes), p(H±) increases with increasing H±, reaches a 
maximum around 5 T, and decreases with further in- 
creasing H± . We define the low (high) field region as the 
field region for which H±_ < 5 T (H± > 5 T), based on 
this behavior of p(H±). This p(H±) dependence gives a 
qualitative explanation of the p(9) dependence [shown in 
the two Insets to Fig. 2(a)] for 9 < 9 C , i.e., at angles at 
which the dissipation is dominated by the perpendicular 
field component H± — Hcos9. Namely, at constant H, 
the resistivity should decrease (increase) with increasing 
9, i.e., decreasing H±, for fields lower (higher) than 5 T 
as, indeed, shown by the bottom (top) Inset to Fig. 2(a) 
for 9 <9 C . 

For the second field orientation (H || a-axis), the p(Hu) 
data of the Inset to Fig. 3 show an initial increase with 
increasing Hu followed by a tendency to saturation for Hu 
approaching 14 T. This monotonic p(Hu) dependence ex- 
plains qualitatively the increase in p(9) [shown in the two 
Insets to Fig. 2(a)] for 9 > 9 C , i.e., at angles at which 
the dissipation is dominated by the parallel field compo- 
nent H\\ = H sm9, in both low and high field regimes. In 
summary, the p(H±) and p(Hu) dependences of the In- 
set to Fig. 3 explain the non-monotonic p{9) dependence 
at low fields and the monotonic p{9) dependence at high 
fields, as shown by bottom and top Inset, respectively, to 
Fig. 2(a). 

Figure 3 shows that the two resistivity curves of its 
Inset, p(H±) and p(H\\), scale for H < 4 T if p(H\\) is 
replaced by p(-f~ 1 H\\), with 7 = 1.7. The physical mean- 
ing of 7 is the anisotropy of the material. The scaling 
of Fig. 3 implies that the same physics is responsible for 
the two scaling behaviors displayed by Figs. 1 and 2 for 
H < 4 T, i.e., fi(Hcos9) and f 2 {H sin 6^, respectively. 

The scaling relationship of Fig. 3 also allows one to 
determine the critical angle 9 C below (above) which the 
resistivity data scale when plotted vs the perpendicu- 
lar (parallel) component of the applied magnetic field 
(Figs. 2(a) and 2(b)). Recall that Eq. (1) gives: 
fi{Hcos9 c ) = f2(Hsin9 c ). This latter relationship and 
the scaling of Fig. 3 [(fi(H) = f 2 (l^ 1 H) for H < 
4 T)] give fi(Hcos9 c ) = h^HsinOc). Therefore, 
9 C = tan^ 1 ^ = 59.5° with 7 = 1.7, as determined above. 
This value of the critical angle is in excellent agreement 
with the 60° value obtained experimentally, which proves 
the consistency of our analysis of the angular and field 
dependent data. So, the critical angle is given by the 
anisotropy of the material. 

The scaling of the magnetoresistivity curves for H\\c- 
axis (for which the Lorentz force is maximum) and for 
H 1 1/| I a-axis (for which the Lorents force is zero) shown 
in Fig. 3 implies that the spin effect, rather than the 
orbital effect, is responsible for the magnetoresistivity in 
region I. The fact that in this region the resistivity is 
also linear in T, indicates that this linear T dependence 
is as well a result of spin fluctuations. This conclusion 
that spin fluctuations dominate the charge transport in 
the non-Fermi liquid regime is consistent with recent In- 
NQR and Co-NMR experiments^, which have revealed 
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that the magnetic nature in CeCoIns is characterized by 
strong antiferromagnetic spin fluctuations in the vecinity 
of the quantum critical point. 

In the following, we show that all the observed fea- 
tures of ADR of CeCoIns in the non-Fermi liquid regime, 
when a magnetic field is rotated from i?||c— axis to 
H 1 1/| \a— axis with 9 the angle between H and the c— axis, 
can be consistently discribed in terms of d-wave density 
waves in a magnetic field. The unconventional density 
wave UDW or d-wave density wave d-DW is a kind of 
density wave in which the gap function A(fc) vanishes on 
the line nodes. Therefore, the transition from the nor- 
mal state to the UDW is a metal to metal transition 29 . 
The UDW exhibits two characteristics: angular depen- 
dent resistivity and giant Nernst effect. A UDW in the 
low temperature phase of a-(ET) 2 KHg(SCN) 4 and the 
metallic phase of (TMTSF) 2 X with X = PF 6 and Re0 4 
has been identified through AD R/Wi? 1 . 

We assume that the electrical conductivity in the non- 
Fermi liquid regime is given by: 



piH.e)- 1 =J2^n{Sech 2 {E n /2k B T)), 



(2) 



where E n is the energy of all fermionic excitations (holes 
and particles). As shown below, this new formula ap- 
pears to be more appropriate than the one used earlier—. 
In the absence of a magnetic field, we assume that the 
quasiparticle energy spectrum is given by: 



E{k) = VC 2 + A 2 (fc), 



(3) 



where £ = v(ku — kp) — (v'/c) cos(ck z ) and A(fc) 
A sin(20) ; here v and v are Fermi velocities in the ab- 
plane and along the c-axis, respectively, fc|| is the radial 
wave vector within the ab— plane, A is the maximum 
value of the d-DW energy gap A(fc), and tan0 = k y /k x . 
Here we assume that A(/c) has d xy -wa.ve symmetry. In 
the vicinity of the nodal points, it is convenient to re- 
place A 2 sin 2 (2</>) by v 2 k\, where k± is perpendicular to 
fen within the ab plane and v 2 /v — A/Ep. Then in a 
magnetic field which makes an angle 9 with the c— axis, 
the energy spectrum becomes^: 



E tn 



±y / 2env 2 H (v\ cos 9\ + v'sinO) — p, 



E± n = ±^2env 2 H (v\ cos 6»| 



p. 



(4) 



(5) 



Here Ef n and Ef n arc the two branches of the Landau 
levels, n = 0, 1, 2,. . ., and fi is the chemical potential. We 
note that in RefiSli only Ef n is considered. The electrical 
conductivity is then given by: 



p(H,t 



ai(Sech 2 (x\ 1 ) + S ech z {x 21 )) , 



»=+,- 



(6) 

where a;^ = Ef 1 /2k B T and xf x = E 2 \/2k B T. For sim- 
plicity, we take only the n = and n = 1 Landau levels. 



Equation (6) fits very well the p(ff) data over the whole 
range of the applied magnetic field, both in the low and 
high field regimes. Representative plots of the resistivity 
data measured in 4, 8, and 10 T along with the fitting 
curves are shown in Fig. 4. In these fits we assume that 



0~n — A n 



1 + I — cos ( 



(7) 



with n — 0, 1. The fitting parameters determined by 
fitting the data measured in different magnetic fields are 
almost constant, hence field independent, as they should 
be. Their values are: the anisotropy 7 = v/v' ~ 1.9, 
y/vv2~ m 1.3 x 10 4 m/s, p = 8.4 K and H = 20 T. This 
value of the anisotropy 7 of 1.9 is very close to the value 
of 1.7 required to scale the data shown in Fig. 3. With 
v = 3.3 x 10 4 m/s, obtained by fitting the H c2 (T) curve 
of CeCoIngiS and the above value of ^vv 2 , we get v 2 /v 
= A(0)/£> w 0.1. Provided that E F = 450 K, this gives 
A(0) ss 45 K, a value expected from the weak-coupling 
theory for d-DW with T c = 20 Kj2S This overall consis- 
tency between the obtained and expected values of the 
fitting paramenters further attests to the appropriateness 
of the model used to analyze the experimental data. 

Moreover, the d-DW model given by Eq. (6) also gives 
the scaling observed experientally at H < 5 T. Figure 5 
is a plot of the p(9) data (open symbols), the fitting with 
Eq. (6) (solid line), and the two scalings (dashed lines) 
calculated from Eq. (6) in which Ef n and E^ n , Eq. (4) 
and (5), respectively, include only the Hcos9 or Hsin9 
terms, for 9 < 60° and 9 > 60°, respectively. Note the 
excellent agreement between the fitting line and the two 
scaling lines. Therefore, even though the model proposed 
here is quite complex, it gives the two experimentally 
observed scaling laws at low magnetic field values. 

In summary, comprehensive angular dependent resis- 
tivity ADR measurements in the non-Fermi liquid region 
of CeCoIns at T < 20 K have shown the presence of 
two different scaling behaviors in the low-field region (re- 
gion I) which involve the H and 9 dependence and are 
due to spin fluctuations. The two scaling regions are 
separated by a critical angle 9 C , which is given by the 
intrinsic anisotropy. In the scaling region, the resistivity 
is linear in T. At higher fields, the ADR data are gov- 
erned by the same physics, even though the two scalings 
fail. A possible explanation for this anomalous ADR in 
terms of c?-wave density waves with Landau quantization 
of the quasiparticle spectrum was presented. This ap- 
proach describes very well the salient features of ADR 
data: (i) the distinctive angular dependences observed 
in the low (dp/dH^ > 0) and high (dp/dHj_ < 0) field 
regimes follow naturally from Eq. (7); (ii) the low field 
scaling behaviors follow from the quasiparticle spectrum 
given by Eqs. (4) and (5). In addition, the quasi two 
dimensional aspect of d-wave density waves is explored 
here for the first time. 
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FIGURE CAPTION 
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FIG. 1: A contour plot of the in-plane resistivity p of CeCoIns 
in the Hcos8 — HsinO plane. Three different regions are 
found in the contour plot, which are distinguished by different 
values of the exponent n of the power-law T dependence of 
the resistivity; i.e., p = po + AT n . The i-axis of this plot is 
along the c— axis of CeCoIns and 8 C is the angle between the 
diagonal of the rectangle in region I and the x-axis (c— axis). 
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FIG. 2: Figure 2(a) Scaling of resistivity p vs the compo- 
nent of the applied magnetic field H perpendicular to the 
afo-planes, Hcos#, measured at a temperature T of 6 K and 
for 1 T < H < 5 T. The arrows mark the points at which 
the curves start to deviate from each other. Insets: p vs the 
angle 9 measured at in an applied mahnetic field of 2 T and 
9 T. 2(b) Resistivity p vs the component of the applied mag- 
netic field H parallel to the afe-planes, Hsin#, measured at 
a temperature T of 6 K in scanning H at fixed orientation 
50° < 9 < 90°. The p vs #sin6> data scale for 60° < 6 < 90°. 
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FIG. 3: Figure 3. In region I, scaling of resistivities p(H±) 
and p(7 _1 ff || ) when plotted vs the component of the magnetic 
field perpendicular to the ab-planes, H±, and parallel to the 
a&-planes, Hu, times the inverse anisotropy 7 _1 , respectively. 
Inset: Plots of p(H±) and p(Hu) curves measured in magnetic 
fields up to 14 T. 
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FIG. 4: Figure 4. Angular (8) dependent resistivity p mea- 
sured at 6 K in a magnetic field of 4 T, 8 T and 10 T. The 
solid lines are fits of Eq. (6) to the data. 
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FIG. 5: Figure 5. Angular (G) dependence of the resistivity p 
measured at 6 K in a magnetic field H = 2 T (open symbols) . 
The solid line is the fit of Eq. (6) to the data. The dashed 
lines are the HcosO and HsmO scalings laws. 



